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Abstract 

Near a solid boundary, E. coli swims in clockwise circular motion. We provide a 
hydrodynamic model for this behavior. We show that circular trajectories are natural 
consequences of force-free and torque-free swimming, and the hydrodynamic interactions 
with the boundary, which also leads to a hydrodynamic trapping of the cells close to 
the surface. We compare the results of the model with experimental data and obtain 
reasonable agreement. In particular, we show that the radius of curvature of the trajectory 
increases with the length of the bacterium body. 



1 Introduction 

The bacterium Escherichia coli (E. coli) has been a micro-organism of choice for studying 
biological and biomechanical processes. In particular, E. coli has been used as the prototypical 
micro-swimmer [1, 2]. 

E. coli and other peritrichously flagellated bacteria swim by the action of rotary motors 
(two to six) embedded in the cell wall. All the motors rotate counter-clockwise or clockwise, 
when viewed from outside the cell, with each motor driving a long, thin, left-handed helical 
flagellar filament. If all the motors rotate counter-clockwise in a viscous fluid (e.g. water), the 
flagella bundle together and propel the bacterial cell forward approximately in a straight line. 
This motion is called a "run". If one or more motors runs clockwise, the flagella unbundle, 
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and the bacteria tumbles. The forward thrust generated by the flagellar bundle in a run is 
opposed by the translational viscous drag on the entire cell. Each flagellum (average length, 
5 to 7 /im) rotates at speeds of approximately 100 Hz [3, 4] and its counter-clockwise rotation 
exerts a net torque on the cell body (average length, 2 to 5 /im). To balance this torque, the 
cell body counter-rotates in a clockwise direction (viewed from behind the organism) at speeds 
of approximately 10 Hz [5]. 

When close to a solid surface, E. coli does not swim in a straight line, but traces out a 
clockwise (when viewed from above the surface), circular trajectory [6, 7, 8, 9, 10, 11]. The 
purpose of this paper is to provide a hydrodynamic model for such circular motion, examine 
how it depends on the size of the cell and other cell parameters, and compare the results of the 
model to experimental data. 

An early observation of circular motion (1971), reported in [7], measured a radius of curva- 
ture for the circles on the order of 25 jum. The swimming direction was clockwise when viewed 
from above, which the authors expected, as the flagellar bundle rotates counter-clockwise and 
the cell body rotates clockwise. The influence of temperature on the motility of E. coli was 
considered in [6]; this work reported circular curves for the motion near a glass slide, with a 
radius on the order of 10 — 50 /im and increasing with temperature. A tracking microscope 
was used later [8, 9] to follow the trajectory of E. coli near a glass surface. Again, near solid 
boundaries, the bacteria were observed to swim in circles, with radius of about 13 /im; the 
authors also found that the swimming speed increased with the distance from the boundary. 
The question of attraction between the swimming bacteria and solid surface has been studied 
in [10, 11], and the distance that cells swim parallel to the surface has been measured (tens of 
nanometers). It was found that standard DLVO (Derjaguin-Landau-Verwey-Overbeek) theory 
could not explain the tendency of the cells to stay near the surfaces, but that some other force 
was still to be identified. The authors proposed that, because of their non-spherical shape, the 
cells swim at an angle to the surface, and therefore constantly swim into the surface. 

Numerically, there has been only one study that has considered the hydrodynamics of a 
swimming bacteria near a no-slip surface [12] (see also early work on flagellar motion near 
boundaries in [13, 14, 15, 16, 17]). The bacterium was modeled as a body of spherical shape with 
a single, solid, helical flagellum and the boundary integral method was used for the numerical 
investigation. In this approach, the total flow field is given by a distribution of fundamental 
singularities for Stokes flow along the surface of the micro-organism. In the simulations, circular 
motion was obtained with a radius of curvature on the order of the length of the micro-organism 
(10 /xm), with a tendency for the micro-organism to swim towards the wall and "crash" into 
it. Furthermore, the authors proposed a physical picture for a clockwise motion. However, no 
simple analytical model was proposed and a numerical integration was required to obtain the 
cell trajectories. 

The goal of this paper is to provide a hydrodynamic model for the motion of E. coli near solid 
boundaries. In §2, we first summarize our experiments to obtain a new set of data on swimming 
speed and circular trajectories for E. coli strain HCB437 near solid surfaces. We then present in 
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Figure 1: Superimposed phase-contrast, video microscopy images show E.coli cells (HCB437) 
swimming in circular trajectories near a glass surface. Left: Superposition of eight seconds 
(240 frames) of video images. Right: Typical superposition of two seconds (60 frames) of video 
images that was used to analyze the length and width of cells, the swimming speed of cells, 
and the radius of curvature of the trajectories. 

§3 our geometrical model for E. coli, and the physical picture for the circular trajectory of the 
bacterium near a no-slip surface, based on the change in hydrodynamic resistance of elements 
along the cell body due to the nearby surface. Using resistive-force theory, we calculate in §4 the 
trajectory of the bacterium. Since the full model requires a matrix inversion to be evaluated, we 
also present an approximate analytical solution for the trajectory. In particular, we show that 
the circular motion is clockwise and that the cells need to swim into the surface as a natural 
consequence of force-free and torque-free swimming. We illustrate the results of our two models 
(the full model and its analytical approximation), show their dependence on various geometrical 
parameters of the cell, and compare the models with our experiments in §5. We find that our 
models are consistent with experimental swimming speeds and radii of curvature of the circular 
motions, and that they allow us to obtain an estimate for the relation between the size of the 
bacterium and its distance to the surface. The values of the various hydrodynamic mobilities 
used in the model are presented in Appendix A, and the cell trajectory far from a surface is 
given in Appendix B. 

2 Experiment 

We examined a dilute suspension of smooth-swimming (i.e. non-tumbling) E. coli cells (HCB437) 
[18] in an observation chamber that was constructed from two glass cover-slips (separation, 80 
/iin) separated by double-sided tape. Videos of cells swimming in circles below the top glass 
surface were collected, digitized, and analyzed. The following parameters were tabulated for 
90 individual cells and their trajectories over a two-second interval: Cell length and width, 
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swimming speed, and the radius of curvature of the trajectory. 

The cells were observed from outside the chamber above the surface, swimming with counter- 
clockwise trajectories; consequently, when viewed from within the liquid (what we will refer to as 
"above the surface" in the remainder of the paper), they are performing clockwise trajectories. 
In Fig. 1 we provide superimposed video images showing the curved trajectories that cells follow 
when swimming near the glass surface. 

2.1 Materials and methods 

2.1.1 Preparation of Motile Cells. 

E. coli strain HCB437 [18] used in these studies is a smooth-swimming strain that is deleted 
for most chemotaxis genes. During cell growth, cells double their length and then divide at 
their approximate midpoint (septate), while maintaining a constant width. The length of 
cells naturally vary depending on the progress of cells through the growth cycle [19]. Media 
components were purchased from Difco or Sigma. Saturated E. coli cultures were grown for 
16 hours in tryptone broth (1% tryptone and 0.5% NaCl) using a rotary shaker (200 rpm) at 
33 °C. Saturated cultures were frozen at —70 °C in 15% glycerol. Motile E. coli cultures were 
obtained by diluting 50 fiL of the thawed saturated culture into 5 mL of fresh tryptone broth, 
and grown in 14 mL sterile, polypropylene tubes at 33 °C on a rotary shaker (150 rpm) for 3.5 
hours. Cells were washed by three successive centrifugations at 2000G for 8 minutes and were 
resuspended into motility buffer [20] (1 mM potassium phosphate, pH 7.0, 0.1 mM Na-EDTA) 
containing 10 mM glucose and 0.18% (w/v) methylcellulose (Methocel 90). Glucose was added 
to maintain motility in an anaerobic environment and methylcellulose was added to reduce the 
tendency of cells to wobble [21]. Filamentous cells were obtained by growing motile cells for 3.5 
hours as described above, adding 50 /ig/mL cephalexin to the culture, and then growing cells 
an additional 0.5 hours [22]. Filamentous cells were then washed as described above. 

2.1.2 Observation of Swimming Cells. 

A volume of 50 fiL of the washed cell suspension (approximately 10 7 cells/mL) was added to 
an observation chamber constructed from two glass coverslips and double-sided tape (Scotch, 
Permanent). The chamber dimensions were approximately: 1 cm wide, 2 cm long, and 80 jum 
high. The microscope coverslips were alternately rinsed with soap and DI water, DI water, 
ethanol, DI water, and then treated with an air plasma for 1 minute at 1-2 Torr (SPI Plasma 
Prep II, power ^80%). The observation chamber was heated to 32 °C using a heated microscope 
stage (Research Instruments Limited). Cells swimming near the upper glass coverslip were 
observed using a Nikon Eclipse E400 upright, phase-contrast microscope. Video images were 
acquired using a 20 x or 40 x Nikon phase objective and a monochrome CCD camera (Marshall 
Electronics V1070) connected to a digital video recorder (Sony GV-D1000) that collected 640 
pixel x 480 pixel images at 30 frames per second. 
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Figure 2: Results of our experimental investigation of swimming E. coli near solid boundaries. 
Left: Radius of curvature of the circular trajectory, 71, as a function of the equivalent sphere 
radius, a, of the elliptical cell body (see text); Right: Swimming speed, U, versus equivalent 
body radius, a. In both cases, we have added the best least-square fit to the data (dashed line). 



2.1.3 Image Analysis. 

Video was captured into a computer using Adobe Premiere and analyzed using ImageJ (avail- 
able for download at http://rsbweb.nih.gov/ij/) or Scion Image (available for download at 
http://www.scioncorp.com) using standard analysis tools. Video images were thresholded so 
that cells appeared black and the background appeared white. The following parameters were 
measured for individual cells in 60 consecutive video frames (2 seconds): The projected area 
of the cell, the midpoint of the cell, and the short and long axis of the cell (approximating the 
cell shape as an ellipse). The average of these values measured over the 2 s interval was used. 
The average cell speed was calculated by measuring the average distance that the midpoint of 
the cell traveled between each video frame and dividing this distance by the video collection 
rate (30 fps or 0.033 s). The radius of curvature of the cell trajectory was calculated by making 
a least square fit of a circle to the two-second trajectory of the midpoint of the cell. A small 
amount of error was introduced by the collection and analysis of cells from multiple regions 
of the swimming chambers and from multiple chambers. Small changes in focus in different 
regions led to variability in the thresholding, which led to some error in the measurement of 
cell widths and lengths. The variability in the cell widths was 0.98 ± .09 /mi; therefore, the 
error due to these small focus changes was less than 10%. 

2.2 Results 

In Fig. 2 we plot the experimental results for the cell swimming speed (U) and the radius of 
curvature of the circles (1Z) as a function of the equivalent sphere radius, a, that is, the radius 
of the sphere that has the same viscous resistance as the prolate ellipsoid of measured cell 
dimensions translating along its axis of symmetry [23]. The large scatter in the experimental 
data, evident in Figs. 1 and 2, can be explained by the natural cell-to-cell variability in the 
number of flagella [24], flagellar rotation rates [4], and distances of cells from the surface [11] 
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Figure 3: Set-up and notations for the mechanical model of E. coli swimming near a solid 
surface. 

(parameters that might also be function of cell length). Nonetheless, the experimental data 
demonstrate that the radius of curvature of the trajectories of cells swimming near a glass 
surface increases with the length of the cell body. 



3 Model 

We present in this section our hydrodynamic model for the motion of E. coli near a flat no-slip 
surface and give a simple physical picture for the circular trajectory. 

3.1 Set-up 

We model the bacterium as a single, left-handed rigid helix attached to a spherical body of 
radius a whose center of mass is located at a distance d above a solid surface, as illustrated 
in Fig. 3; the liquid gap between the solid surface and the cell body has height h. The cell 
is assumed to be parallel to the surface and oriented in the ^/-direction. The helix is assumed 
to have thickness 2r, radius b, wavelength A, with a number n of wavelengths along the helix, 
such that the total length of the helix along the ^/-direction is Lii = nX. The assumption of 
sphericity, although not completely realistic for the cell body of E. coli, which is more like a 
2:1 prolate ellipsoid, was made in order to use well-known mobility formulae, and we expect 
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Figure 4: Physical picture for the out-of-plane rotation of the bacterium: (a) The positive 
y- rotation of the cell body leads to a net viscous x- force on the cell body, T], > 0; (b) The 
negative y-rotation of the helical bundle leads to a net negative viscous x-force on the nagella, 
T^. < 0. The spatial distribution of these forces leads to a negative z-torque on the bacterium, 
which makes it rotate clockwise around the z-axis. Therefore, when viewed from above, the 
bacterium swims to its right. 

therefore our results to be correct within a shape factor of order unity. Due to the action of 
rotary motors, the bundle is rotating in the counter-clockwise direction (viewed from behind) 
with an angular velocity ui = —uj e y relative to the body, with uj > (see Fig. 3). We denote by 
U = (U x , Uy, U z ) and ft = (Q x , fl y , Q z ) the instantaneous velocity and rotation rate (measured 
from the center of the cell body), respectively, of the bacterium. 

3.2 Physical picture 

In the absence of a nearby wall, the bacterium swims approximately in a straight line, U = U y e y , 
and rotates along its swimming axis, O = Q y e y . The velocity U y > is obtained by balancing 
the propulsive force of the helical bundle with the viscous resistance on the whole bacterium 
and the rotation rate VL y > is found by the balance of viscous moments about the y-axis (see 
Appendix B). 

What changes when the micro-organism is swimming near a solid surface? The cell body 
and the helical bundle contribute together to a rotation of the bacterium around the z-axis 
(see notation in Fig. 3) (see also [25]). First, as the cell body is near the surface, when it 
rotates about the y-axis at a rate Q y > 0, there is a viscous force acting on the cell body in the 
x-direction, T^e x , with T]. > (this is a standard hydrodynamic result, see [26]). Furthermore, 
the bundle of nagella is also acted upon by a net force in the x-direction. Since the bundle takes 
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the shape of a helix, parts of the bundle are located close to the surface and parts are located 
further away. The local drag coefficient on an elongated filament decreases with increasing 
distance to the nearby surface (see details below), which means that the parts of the bundle 
that are close to the surface will be subjected to a larger viscous force compared to portions 
of the helix located further away from the surface. As the bundle rotates counter-clockwise 
around the y-axis (viewed from behind), the portions of the helix that are closer to the surface 
have a positive x-velocity, and therefore the net viscous force acting on the bundle, J^e x , is 
ncg ative, Tl < 0. Of course, since the swimming bacterium as a whole is force- free, we have 
necessarily = —T].. 

As a consequence of the viscous forces acting on both the helical bundle and the cell body and 
their spatial distribution, a negative torque, C z < 0, will act on the bacterium and will rotate 
the entire cell clockwise around the z-axis (Fig. 4). When viewed from above, the bacterium 
will therefore swim to the right, as is observed experimentally. Note that, since the bacterium 
as a whole is torque-free (the Reynolds number is low, Re ~ 10~ 4 , so forces and torques need 
to balance at each instant), this torque will be balanced by a positive torque arising from the 
viscous resistance to a rotation around the z-axis. 

This physical picture allows us to obtain an estimate for the radius of curvature 1Z of the 
motion, as the ratio of the swimming velocity U y to the out-of-plane rotation rate Q z . Since 
the Reynolds number is low, the equations of motion for the fluid are linear (Stokes flow) and 
therefore instantaneous viscous forces and torques for various parts of the bacterium are linearly 
related to their velocities and rotation rates. Let us denote by M. and M the viscous (tensorial) 
mobilities of the bacterium flagella and body respectively which are non-zero away from a solid 
surface, and denote by W and V those mobilities that arise due to the presence of the nearby 
surface and tend to zero far from a wall (see §4 for details). For all these mobilities, we will 
use notations of the form Ai^f, where the superscript a/3 is either FU (in which case M-f^ 
denotes how the zth component of a viscous force is linearly related to the jth component of 
a velocity), FQ (force - rotation rate), LU (torque - velocity) or LQ (torque - rotation rate). 
We also always use the convention that the mobilities are positive, and might therefore appear 
with a minus sign when necessary. 

The swimming velocity is obtained by balancing the propulsive force of the micro-organism, 
A4™(u — Q y ), with the viscous drag on the bacteria (M-yy + ■^ y F y u )U y , so that 



The rotation rate can be estimated by balancing the wall-induced torque mentioned above, 
C z ~ W^iuj — fly), with the viscous resistance to rotation of the whole bacterium, which is 
mostly due to the viscous resistance of the flagella, — M^fl z , that is 



The two previous balances lead to an estimate for the radius of the circular motion as 



M™(u - n y ) « (M 




(1) 



(2) 



Uy M™M 



(3) 



s 



As is demonstrated below, this simple estimate is consistent with a more detailed calculation, 
as well as with experimental results. 



4 Trajectory calculation for the bacterium 

We proceed in this section by presenting the detailed calculation for the trajectory of the 
bacterium using resistive-force theory for the flagellar hydrodynamics, and exploit it to obtain 
an approximate analytical solution. 



4.1 Modeling of flagella hydrodynamics 

The modeling chosen here for the helical hydrodynamics is that of resistive- force theory (RFT), 
as first introduced by Gray and Hancock [27], since it is the simplest approach to the zero- 
Reynolds-number hydrodynamics of elongated bodies. The method is an approximation to the 
equations of slender-body theory (SBT). SBT considers the zero- Reynolds- number dynamics 
of long and slender filaments by distributing fundamental Stokes flow singularities at their 
centerline (Stokeslet and force-dipole) [26, 28]. The idea was first introduced by Hancock [29], 
is reviewed in detail by Lighthill [30], and has been applied to the case of helical flagella in [31]. 

RFT is the leading-order approximation of SBT, which gives results accurate at order 
(9([log(L/r)] _1 ), where L is the length along the filament and r its radius; the next term 
in the formulation of SBT appears at order 0(1). The complexity of fully solving for the spa- 
tial distribution of singularities on a moving flagellar filament is replaced by introducing a set 
of local drag coefficients. Let us consider a portion of the filament of length 8£, oriented along 
the tangential vector, t, and moving at a velocity u in a viscous liquid. The local velocity can 
be decomposed into a parallel and perpendicular components, u = uy +ui, where uy is parallel 
to the tangential vector, uy = (u-t)t, and uj_ is perpendicular to it, uj_ = u — U| . RFT assigns 
values for the local drag coefficients, cm and c±, which relate the local viscous force per unit 
length to the local parallel and perpendicular velocities, such that the total force on an element 
of length 5£ can be written 

<SF = -^(c[|U|| + c x u.i). (4) 

For a periodic flagellar filament (wavelength A) performing planar oscillations far from a solid 
surface, we have approximately [27, 17] 

Cj_ = 2 Cm- (5) 



11 ln(2A/r) - 1/2 : 

The case of helical flagella was first considered in this context in [32]. Note that the drag 
anisotropy between tangential and perpendicular motion is the fundamental origin of the flag- 
ellar propulsion of micro-organisms [2, 27, 30]. Although it is only an approximate method, 
RFT has been shown in the past to provide both qualitative and quantitative information about 
the locomotion of micro-organisms [17, 27, 30, 33, 34]. 
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The presence of a solid surface modifies the values of the resistance coefficients for both the 
cell body and its flagella [35, 36, 37, 38, 39, 17, 15, 16, 14]. Elements of the helical flagella are 
located at a distance d(z) ranging between d — b and d + b to the solid surface, which are both 
smaller than the helix wavelength A, so that the viscous resistance to motion of the flagella 
is dominated by the interactions with the surface. Since r <C d ± b, we consider the far-field 
asymptotic results of [15] (see also the review in [17]) and use 

27T?7 



ln(2d(z)/r)' 



c± = 2 cm 



(6) 



Deviations from 2 for the ratio c±/c» were discussed in this context by Katz & Blake [16]. We 
will denote by cj[ the value of the drag coefficient, Eq. (6), when d(z) = d, and will denote 
deviations from this value by the function /, so that c\\(z) = c\\f{z). 

4.2 Mobilities 

We consider separately the mobilities of the cell body and its flagella, neglecting therefore the 
hydrodynamic interactions between these two parts of the micro-organism. Although this an 
approximation, we expect it will contribute only to a small error in the final results as the 
presence of a nearby surface leads to spatially localized flow fields, decaying at least as fast as 
a Stokeslet-dipole (~ 1/r 2 ). 

As described earlier, we denote by M. and M the mobilities that are non-zero even in the 
absence of a wall, and by V and W those which are equal to zero when the micro-organism 
swims far from the surface (with the conventions that the mobilities are positive). The mobility 
matrix for the spherical cell can be written as 
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with values calculated in Appendix A. As can be seen in Eq. (8), the matrix B is almost full; 
the elements reported to be zero are either exactly zero at all instants or time-average to zero 
over the rotation period of a flagellar filament, T = 2-k/uj. If we define 

x = {u x ,u y ,u g ,n x ,n v ,n z ) T , and y = (0,0,0, o,^,o) T , (9) 

then the requirement that the micro-organism is free-swimming, C = and J- = 0, becomes a 
6x6 linear system to solve for X of the form 

(A + B)X = By. (10) 

The solution for the velocity, U, and the rotation rate, fl, can be found by simply substituting 
the values of the mobilities from Appendix A and solving the linear system Eq. (10). The radius 
of curvature of the in-plane motion will then be given by 

U 



a 



n = —, u = jm + m- (111 



V 



4.3 Approximate analytical solution 

When the bacteria swims far from the solid surface, an analytical solution to motion can be 
found and we give it in Appendix B. In the presence of a solid surface, an analytical solution to 
the linear system, Eq. (10), exists in theory by direct matrix inversion, but it is very complicated 
and not very enlightening. We present below instead an approximate analytical solution of the 
linear system. 

First, we note that, in the case of E. coli, a number mobilities can be neglected between the 
elements of A and B. They are 

AC«AC, M£<:Mf£, (12a) 

M™ « M™, M% « A^ n , (12b) 
Wxy « V%, « V™, (12c) 

W™ « V™ W™ « V™. (12d) 

Furthermore, since the x- and ^-components of both velocity and rotation rate are zero 
far from the solid surface, we make the assumption that, near the surface, these components 
are at most on the order of the j/-components: we therefore assume that (U x , U z ) < U y and 
(Q X ,Q Z ) < Q y . We further assume that Q y u, as is the case far from the surface. Finally, 
since we have in general (W°f , V"f) (-A^ J\f? ), where j = x or z, these assumptions allow 
us to simplify further the mobilities in the matrices A and B. 

In that case, the equations C y = and Ty = lead to the approximate solutions for 
the swimming speed and body rotation 

M FQ 

U V « M FU fjfFU ( 13& ) 

J l yy ^-IVyy 

n „ KIKl u (13h) 
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and Qy is indeed verified to be much smaller than u. We can then use ^2 T z = and obtain 



(14) 



It follows, by substituting Eq. (14) into C z = and evaluating the leading-order contribution, 
that 

1 



U.r 



M LU M FQ 
J l zz J l zx 
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* v y y 



q x - M^n 



(15) 



As a consequence, substituting Eqs. (14) and (15) into C x = 0, using Eq. (13a) and evaluating 
the leading-order term leads to 



M™n x + — f , TTT n z + 
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Finally, substituting Eqs. (13a), (13b) and (15) into Y^-^x = and keeping the leading-order 
terms leads to 
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Solving the 2x2 linear systems of equations given by Eqs. (16)-(17), and keeping only the 
leading-order terms leads to approximate formulae for the x- and ^-components of the rotation 
rates as 



n, 



V LU M FQ 

xy J yy 



zz M F ¥ + N FU 



UJ. 



(18a) 
18b) 



Note that the denominator in the equation for fl z , Eq. (18b), is dominated by M.^ but not by 
much, so we need to keep both terms to obtain correct orders of magnitude. These equations 
allow us to verify that, for E. coli, Q x is much smaller than Q y and Q z is of the same order as 
Q y . Note also that we obtain Q z < 0, which means that the bacteria is swimming to its right 
(clockwise trajectory viewed from above) and that Q x < 0, so that the bacteria will also have 
the tendency to swim into the surface. Also, we observe that 



<1, 



(19) 



so the time scale for re-orientation of the bacteria perpendicular to the surface is much larger 
than the typical swimming time scale. 
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Table 1: Summary of the results of the simplified model for E. coli swimming near a solid 
surface. The mobilities are calculated in Appendix A. 

Now, substituting Eq. (18a) and Eq. (18b) into Eq. (14) and Eq. (15) and keeping leading- 
order terms leads to 

v lu M Fn m fci 

tt ~ x v zx yy (onai 

y^LQ 

Ux * ^aoCT) T^ Ut (20b) 



- Mi 



and we get that U x > and, more important, that U z < 0. This result, together with the result 
that Q x < 0, shows that hydrodynamic interactions "trap" the cell close to the wall. Note that 
this trapping does not require cells to be non-spherical. Note also that 

yy V zx J 

where e = 2nb/X and J is defined in Appendix A, and 



so the calculation assumptions are consistent. 
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Figure 5: Comparison between the results of the experiments (>), the full hydrodynamic model 
(Eq. 10, ■. and best fit, straight line) and the simplified model (Table 1, dashed line). Top: 
L\\ = 5 fim, and the values leading to the best least-square fits between the experiments and 
the full hydrodynamic model, h = 38 nm and ou — 211 Hz: (a) Radius of curvature, 1Z, and 
(b) swimming velocity, U, as a function of the bacterial radius a. Bottom: Same as top but for 
L\\ = 7 /im, h — 10 nm and uj = 194 Hz: (c) Radius of curvature, and (d) swimming velocity 
as a function of the bacteria radius a. 

We can finally evaluate the approximate solution for the radius of curvature of the circular 
trajectory It is given by 

v = JL „ 3l „ M™M™ ( _ M™M% 

\n z \ ~ |n,| ~ y^(Mff + M y F y u ) V m™(m™ + N™) 

which is very similar to that given by the simple physical picture, Eq. (3). 

The results of the analytical model are summarized in Table 1. When we set V = W = 0, 
and assume that the previous approximations still hold, the results from Appendix B (swimming 
far from surface) are recovered. 
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5 Results of the model and comparison with experi- 
ments 



We present in this section the comparison between the results of our hydrodynamic model 
and our experiments. The geometric characteristics of the flagellar bundles that we use are 
r = 20 nm, A = 2.5 /mi and b = 250 nm [5, 4]. The length of the flagella is usually in the range 
from L|| = 5 to 7 /mi, and we will test both values. For the cell radius a, we take the equivalent 
sphere radius a that has the same viscous resistance as the prolate ellipsoid of measured cell 
dimensions translating along its axis of symmetry [23]; we use however the measured value of 
the minor axis of the elliptical head to correctly estimate the distance between the attached 
helical flagella and the wall and this results in a small scatter in the theoretical predictions (■) 
of Fig. 5. The only parameter in the model whose value is unknown is the gap thickness h. The 
minimum distance cells can swim from the surface is about 10 nm because of the protrusion of 
the flagellar hook from the cell body. Values of h have been measured to be 30 — 40 nm [1 1] . 
In order to compare the model with our experimental data, we will assume h to be in the range 
from 10 to 50 nm. 

Despite the large scatter in our experimental data, we find that the results of our hydro- 
dynamic model are consistent with our data, both for the value of the radius of curvature of 
the trajectory, 1Z, and the swimming speed, U = {U 2 X + U^ 2 . Typical results comparing 
experiment and theory are illustrated in Fig. 5 for Ly =5 /mi and 7 /mi. The values of the 
distance to the wall h and the flagella rotation speed io were chosen to lead to the best fit of 
radius of curvature and velocity of the full hydrodynamic model (square symbols and straight 
line) to the experimental data (triangular symbols); we find h = 38 nm and to = 211 Hz for 
Lii — 5 /mi, and h — 10 nm and u = 194 Hz for Ln = 7 /mi. These values are consistent with 
the measurements of [11] and with typical values for the rotation rate of flagella in E. coli [3]. 

Furthermore, we observe in Fig. 5 that the approximate analytical model (Table 1, dashed 
line) gives a good approximation of the swimming speed for both values of L\\ but fails to predict 
the correct radius of curvature of the trajectory for L\\ = 7 /im (the dashed line appears however 
to go to through some of the experimental data due to the large scatter in the measured radius 
of curvature). More importantly, the approximate analytical model fails to capture the increase 
of the radius of curvature of the trajectory with the cell size, an experimental feature that is 
predicted by the full hydrodynamic model. 

Let us now investigate the dependence of the models, that is the values of {7Z, U } given by 
both the full model and the approximate analytical model, with other geometrical parameters 
describing the bacteria: {b, h, r, A} 1 . In order to display the variations, we will fix the values 
to be r = 20 nm, A = 2.5 /im, b = 250 nm, Lu = 7 /im, h = 30 nm, a = 1 /mi and u = 200 Hz, 
and will then vary one of these parameters at a time. The results are displayed in Fig. 6 for the 
full hydrodynamic model (solid lines, squares) and the approximate analytical model (dashed 

1 The dependence of the model on the cell body, a, is illustrated in Fig. 5. Moreover, from Eq. (10), it is 
straightforward to see that both U and f2 scale with cu, and therefore TZ is independent of u. 
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Figure 6: Dependence of the results {U, 1Z} on the geometrical parameters {b, h,r, A} for the 
two models (full model: squares and straight line; approximate analytical model: circles and 
dashed line), in the case where r = 20 nm, A = 2.5 /im, b = 250 nm, L\\ = 7 /iin, h = 30 nm, 
a = 1 /im and ui = 200 Hz, and one of the parameters is varied at a time. Dependence on 
(a) the helix radius b, (b) the gap thickness h, (c) the helix half-thickness r, and (d) the helix 
wavelength, A. 



line, circles). 

These results first confirm that both models are in agreement for the trends and values of 
the swimming velocity, but that the approximate analytical model overestimates the value of 
the radius of curvature of the trajectory (by up to 50%). The dependence of the swimming 
velocity, U, is found to be consistent with the increase of the propulsive viscous force with 
b/X and r, and the decrease of the viscous resistance with h (see the values of the mobilities 
as calculated in Appendix A). The radius of curvature increases with h, consistent with a 
decreases of the induced ^-torque with the distance to the surface. The radius also decreases 
with r, consistent with an increase in the hydrodynamic interactions with the nearby surface 
as described by Eq. (6). Furthermore, 1Z decreases with b, confirming the important role of 
the viscous resistance on parts of the helix that are close to the surface (whose distance to the 
surface decreases with b) in inducing the torque on the cell in the z-direction. 

Finally, we note that the value of the radius of curvature from the model depends strongly on 
the unknown gap thickness h. Returning to the comparison with the results of our experiments, 
we see that data for larger cells tend to be more consistent with the model for large values of 
h (Figs. 5 and 6). Thus we propose here that, if we suppose that all bacteria have the same 



1(3 




a (n m) 

Figure 7: Best fit to the experimental data (>) by a h(a) law for L» = 7 fim and u = 194 Hz, 
as given by Eq. (24) (■ and best fit, straight line). The dotted line represent the best fit of 
the experimental data by a straight line (in the least-square sense), which we have added as a 
guide to the eyes. 

geometrical characteristics, {b,r,\,L\\}, our hydrodynamic model could be used to estimate 
the relation between the typical cell size, a, and its steady-state distance to the wall, h, that 
is, the h(a) relationship. The results are illustrated in Fig. 7 where we have plotted together 
the results of the experiments with the predictions of the hydrodynamic model with a linear 
relationship between a and h, 

Ka) = h + hi, (24) 

which leads to the best overall fit to the experimental data; the parameters for this fit are 
ho = 10 nm, a\ = 0.7 /im, a<i = 0.8 /xm and hi = 11 nm. 

6 Conclusion 

We have presented a hydrodynamic model for the swimming of E. coli near solid boundaries 
and compared it to a new set of measurements of cell velocities and trajectories. We have shown 
that force-free and torque-free swimming was responsible for the clockwise circular motion of 
the cells, as well as for their hydrodynamic "trapping" close to the surface, that is, Q x < and 
U z < 0. The fact that cells do not eventually come in contact with the surface is probably due 
to other short-range surface-cell interactions [10, 11]. 

The main assumptions made in this paper, and which illustrate the differences between 
real swimming E. coli cells and our model, are the following: (1) We have replaced the bundle 
of several flagella by a single rigid helix; (2) We have assumed that the the cell body was 
spherical; (3) We have ignored all interactions between the cell body and the flagella. Relaxing 
these assumptions would improve on the agreement between theory and experiments, but we 
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do not expect it would change the physical picture given in this paper. Including the presence 
of a second (top) boundary should also modify the cell trajectories. 

Finally, we wish to remark that if the surface was a perfectly-slipping interface (such as the 
free surface between air and water) instead of a no-slip surface, the change of the direction of 
the image system for a point force [40] should lead to bacteria swimming in circles, but in a 
counterclockwise direction (X. L. Lu, University of Pittsburgh, private communication). 
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A Cell mobilities 

We present in this Appendix the values of the hydrodynamic mobilities of the bacteria. First, 
since we have h <C a, the lubrication approximation can be made to derive the mobilities for 
the cell body [35, 36, 37, 38]. We find that they are given by 
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Note that we assumed that Af^f 1 was equal to its far-field value, as it was shown that the 
presence of a nearby surface has only a small effect on the value of this mobility [35]. 
Second, the bundle of helical flagella is described by the equation 
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6sin(s — ut) 
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where s ranges from to 2nir J and u is the rotation rate of the flagella bundle relative to the 
cell body. In that case, the mobility calculation was done according to the principle introduced 
in §4.1 and we get 

M™ = Ae = 2^,1+^, (27a) 



(1 + e 2 )V2 



1 + 2e 



2 



Kv = (1 + e2) i/ 2 » ( 27b ) 

Kv = M^^bL ]{ (1+ £ e2)1/2 , (27c) 

= M™ = ^L\±0^ 2 , (27d) 

M™ = M™ = M% = M% = ^L\ , (27e) 



A4f? = 2c\\tfh ji ' (27f) 



(1 + e 2 )V2 



AC = = -M™ = M™ = icjf 6L|, (1 + £ e2)1/2 , (27g) 

VC = Wy = -q|L B 7 - * 1/2 J, (271) 



11 (1 _|_ e 2^1/2" 

e 

(l + e 2 ) 1 ^ 



wJT = vC = -qi L Li , ! 2 m/ 2 ^ (27j) 



= V^P = -2c|^! (i + ^iA ^ (271) 
>V^ = W^ = -q6 2 L|| (1+ £ e2)1/2 X, (27m) 

where e = 2nb/X, and where we have defined the two integrals 

X — / cos(27TM)/(cos(27ru))du, 3= \ (u + u ) cos(27rmt)/(cos(27rmt))du. (28) 
Jo Jo 

Note that for the calculation of M.yy, the contribution due to the local rotation of the flagella 
can be neglected because r <t^b [32]. 
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B Swimming far from a surface 



matrices become 
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Solving Eq. (10) for the velocities and rotation rates in this case leads to U x = U z 
and 
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In the absence of a wall, the bacteria swims therefore in a straight line and rotates its body in 
the direction opposed to that of the flagella. 
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